Modeling noise induced resonance in an excitable system: An alternative approach 
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Recently, it is observed [Md. Nurujjaman et al, Phy. Rev. E 80, 015201 (R) (2009)] that in an 
excitable system, one can maintain noise induced coherency in the coherence resonance by blocking 
the destructive effect of the noise on the system at higher noise level. This phenomenon of constant 
coherence resonance (CCR) cannot be explained by the existing way of simulation of the model 
equations of an excitable system with added noise. In this paper, we have proposed a general model 
which explains the noise induced resonance phenomenon CCR as well as coherence resonance (CR) 
and stochastic resonance (SR). The simulation has been carried out considering the basic mechanism 
of noise induced resonance phenomena: noise only perturbs the system control parameter to excite 
coherent oscillations, taking proper precautions so that the destructive effect of noise does not affect 
the system. In this approach, the CR has been obtained from the interference between the system 
output and noise, and the SR has been obtained by adding noise and a subthreshold signal. This 
also explains the observation of the frequency shift of coherent oscillations in the CCR with noise 
level. 



I. INTRODUCTION 

Noise induced resonances: coherence resonance (CR) 
and stochastic resonance (SR) in a threshold or excitable 
systems, have been studied both numerically and experi- 
mentally in many physical, chemical, biological and elec- 
tronics systems^ B i i, i, i, 0, 1 i 0, (ll E 13 Q 
IitI . ITsI . Il^ . The regularity of the dynamical be- 
havior of an excitable system emerges by virtue of the in- 
terplay between the autonomous nonlinear dynamics and 
the optimum superimposed noise that is termed as CR. 
Whereas, in SR, the response of the system to a weak pe- 
riodic input signal is amplified or optimized by the pres- 
ence of a particular level of noise. Recently, another type 
of CR has been observed in plasma and electronic sys- 
tems, where the coherency remains almost constant even 
at higher noise level [1, 0, Q • This CR phenomenon of 
constant coherency may be termed as constant coherence 
resonance (CCR). The main difference between CR and 
CCR is that, in CR autonomous dynamics is destroyed 
at higher noise level, whereas in CCR, system dynamics 
remains unaffected from the destructive effect of noise. 
One of the experimentally observed feature of the CCR 
is that the frequency of coherent oscillations increases or 
decreases with increase in noise level depending on au- 
tonomous dynamics of the system. So far there is no 
theoretical explanation of the CCR and frequency shift. 

Now the qualitative features of an excitable system 
can be obtained using two nonlinear autonomous dif- 
ferential equations which have been used for modeling 
the noise induced resonances (23| . In the existing liter- 
atures, the CR phenomenon has been modeled just by 
adding noise to any one of the above equations, where 
maximum coherency has appeared at optimum noise 
level [H, mi, mm [2i|. But in this approach, one can- 
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not explain the appearance of CCR. Even explanation 
of SR is not very clear, though there are several experi- 
ments on this phenomenon ^, .26,] . The main problem of 
adding noise directly to the equations is that the noise 
not only perturbs the system control parameter, but also 
destroys the actual dynamics. The effect of adding noise 
directly to the equations has been discussed in the con- 
text of the FitzHugh-Nagumo model in Section IIVI In 
this paper, our goal is to present a general model which 
explains CR, SR and CCR. The model has been devel- 
oped based on the excitation mechanism of the system 
using FitzHugh-Nagumo model, a well known paradigm 
for modeling noise induced resonances. 

Rest of the paper has been organized as follows: we 
have discussed the excitation mechanism of an excitable 
system and effect of noise in Section [TTl Based on 
the physical arguments of Section |TT1 we have formu- 
lated the model for the noise invoked resonances in Sec- 
tion [nTl The results and discussion of the simulation us- 
ing FitzHugh-Nagumo has been presented in Section [TVl 
Finally a conclusion has been drawn in Section [V] 



II. EXCITABLE DYNAMICS AND EFFECT OF 
NOISE 

The basic characteristic of an excitable system is that 
it shows a fixed point or coherent limit cycle oscillations 
depending on the value of the control parameter (CP) 
of the dynamics. The point where the system changes 
from oscillatory to fixed point behavior is called thresh- 
old or bifurcation point. Now if the system shows fixed 
point behavior for the value of CP below the threshold, 
then the dynamics will be limit cycle oscillations on the 
other side of the threshold or vice versa, depending upon 
the system properties. The frequency of the autonomous 
oscillations in the excited state changes with control pa- 
rameter. For example, frequency of the oscillations in 
FitzHugh-Nagumo model increases when one increases 



2 



the control parameter from the threshold value. The in- 
crease in frequency of the coherent oscillations has also 
been observed in real experiments [3, 0, 0- This has 
important effect on the increase in the frequency of the 
coherent oscillations in CCR with increase in noise level, 
that has been discussed at the last of next paragraph. Re- 
sponse of an excitable system at fixed point state to an 
external perturbation applied on the CP depends on the 
perturbation amplitude. When the amplitude is small 
such that the threshold is not crossed, the system remains 
at its fixed point state, and when it is large enough to 
cross the threshold, the system returns to its fixed point 
deterministically, i.e., once the threshold is crossed, the 
system becomes almost independent of the perturbation 
and comes to its fixed point state traversing one limit 
cycle [il,i[iia[2li. 

For stochastic perturbation (i.e., noise), limit cycles 
will appear randomly due to random crossing of the 
threshold. The time between appearance of the consec- 
utive limit cycles can be splitted into two characteristic 
times: an activation time (ta) through which system re- 
sides at the fixed point state between two consecutive 
limit cycles, plus a refractory time (tr) which is taken 
to traverse one limit cycle. Activation time {ta) will 
decrease gradually with noise, as the threshold will be 
crossed more frequently at higher noise amplitude. When 
the noise amplitude is fairly large and its changing fre- 
quency is higher than refractory frequency (which is gen- 
erally the case in experiments), then the threshold will be 
crossed several times during one limit cycle that makes 
ta « and the inter-peak distances are practically deter- 
mined by the time period (t,.) of the limit cycle. At this 
stage, the dynamics is almost similar to the coherent os- 
cillations at the excited state and is termed as coherence 
resonance (CR). From this stage, increase in the noise 
may lead to two possible effects: firstly, it destroys the 
structures of the limit cycles or the actual dynamics of 
the system, that leads to decrease in the coherency, which 
has been observed in many experiments [1, 0, S S 0, H, 
i, [M HH, E M M MMMMM- secondly, it per- 
turbs the CP across the threshold so frequently that the 
system remains at its excited state and produces almost 
coherent oscillations, i.e., constant coherence resonance 
(CCR). CCR is not always realizable in the real experi- 
ments, as at higher noise level it is difficult to keep the 
noise away from its destructive effect. CCR has been 
observed in some recent experiments 0, 0, ■ Another 
interesting experimental observation regarding the CCR 
is that the frequency of the noise invoked oscillations 
increases with increase in noise level in those excitable 
systems whose frequency of the autonomous limit cycle 
oscillations increases with increase in CP from the thresh- 
old iiSi. 

Now we can summarize the effects of noise on excitable 
system: (a) Noise helps the system to cross the threshold 
by perturbing the CP; (b) it may destroy the system 
dynamics at higher noise level; (c) If the destructive effect 
of noise is blocked, the system may remain at excited 



state and hence may show CCR; (d) In case of CCR 
noise is also responsible for the effective increase in CP 
above the threshold. All these facts have been used to 
formulate a general model which has been described in 
the next section. 



III. FORMULATION OF THE MODEL 

The dynamical state of an excitable system can be 
represented by a two dimensional state vector X = 
{x, nY' € and the equation of motion in general form 
can be described by two first order differential equa- 
tions ^^j^ — F{X{t),a), whose fixed point is given by 
^ = = F{X) [1^, where, F is suitable nonlinear func- 
tion and a is the control parameter. Let a = ath be the 
bifurcation point (or threshold) of the system at which it 
goes from the fixed point to the oscillatory state. Now, 
the value of a = ao is so chosen that the system remains 
in the excitable fixed point region. Then the dynamical 
behaviour of the above system with Gaussian noise (^), 
can be written as 

= F(X(i),ao + ^0; for + ao > ath 
= 0; for D^ + ao < ath 

where, represents the amplitude of the noise of 
standard deviation D. Here the simulation has not been 
carried out just adding random noise at each step of in- 
tegration as was done in the previous works. The in- 
tegration has been carried out based on the physics of 
an excitable system discussed in the previous section. 
Whenever CP crosses the threshold due to noise pertur- 
bation, system traverses one limit cycle. Now once the 
threshold is crossed, the noise remains ineffective in the 
dynamics until the system comes to its fixed point state, 
even if the CP crosses the threshold several times during 
this refractory phase, i.e., at the time of traversing one 
limit cycle 123| . This indicates that the noise need not be 
added at every step of integration during the simulation. 
Hence CP should be kept constant up to one limit cycle 
during the integration. This is also clear from the points 
(a) and (c) discussed in the previous section [Section |TT] . 
The possible value of the CP for each limit cycle inte- 
gration can be guessed from the point (d) [Section |TT] . If 
the noise would only lift the system at its excited state, 
then the frequency of the noise invoked dynamics would 
not change with increase in the noise level. But from the 
experimental results, it is clear that the frequency of the 
oscillations increases with increase in the noise level. So 
only possible reason for increase in the frequency is due 
to the increase in the effective value of CP. Hence, one 
can guess that the effective value of CP at the start of 
the limit cycle can be taken to be its value at fixed point 
state plus the noise amplitude at the time of excitation 
(which remains constant during one limit cycle). 

Now for Eqns.[l] whenever, the condition D^+ao > ath 
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is satisfied, the system will be lifted to its excited state 
and DS^ + aQ will remain constant during each limit cycle. 
During this time the system does not feel the effect of 
noise, i.e., it remains independent of noise perturbation. 
When it comes to its rest state, it feels again the noise 
perturbation and whenever, the condition + oq > ath 
is satisfied, the system traverse one limit cycle with a = 
ao + D^. By imposing this condition, noise has been 
blocked from its destructive effect on the system. In some 
cases like FitzHugh-Nagumo model, the system shows 
oscillatory behaviour when a < ath- In this case the 
condition to be satisfied by noise is + oq < ath which 
will be satisfied by the negative amplitudes of the noise. 
But the main mechanism is same for both the cases. 



IV. RESULTS AND DISCUSSION 

There may be several choice of F [Eqn. [T] so that it 
shows an excitable dynamics (23l |. In the present paper, 
the famous two dimensional Fitz Hugh-Nagumo model 
for excitable system has been studied, whose equations 
of motion are [22] 
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FIG. 1: Frequency of the limit cycle oscillations with {ath — a) 
from FFT (solid line) and from inverse of the refractory time 
(dashed line). 

where, e — 0.01 and the control parameter 'a' govern 
the dynamics. For |a| > 1 and \a\ < 1, the system [Eqn [2] 
shows fixed point and limit cycle oscillations respectively. 
Therefore, ath = 1 is the threshold. Hence, the dynamics 
is fixed point and oscillatory above and below this point 
respectively. Frequency of the limit cycle oscillations in- 
creases between < a < 1 and then decreases between 
> a > — 1. In Fig[Tl the solid and dashed lines show 
increase in the frequency obtained using FFT and of 
the oscillations for parameter 1 < a < 0.5 respectively 
and both estimations are identical. It also shows that the 



frequency of the autonomous dynamics increases with in- 
crease in CP below the threshold. In the next subsection 
we will see that the increase in frequency of the coherent 
oscillation in CCR follow the trend due to noise. 



A. Coherence and Constant Coherence Resonance 

To study the noise invoked dynamics, 'a' is set to 
ao = 1.05 so that Eqn [2] shows fixed point behavior and 
is perturbed by noise (^) of standard deviation D. When- 
ever, a — ao + D^ < 1 is satisfied, the system will traverse 
one limit cycle before coming to its fixed point state. The 
refractory time i,- will depend upon the parameter value 
a = ao + as shown in Fig [TJ For low level of noise, 
the limit cycle appears sparsely and increases with in- 
crease in noise level. The regularity of the appearance 
of these oscillations is determined by coherence param- 
eter NV = a(Ti)/mean{Ti), where a is the standard 
deviations of the interpeak distances (T^) of the limit cy- 
cles [13 . Least NV indicates high degree of coherency 
in the system output. 
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FIG. 2: Dashed line shows the NV of the limit cycles when 
the destructive effect of noise is blocked. The solid line shows 
the NV computed after interfering the output with applied 
noise. 

Fig [2] shows NV vs D plot for two different cases. The 
dashed line [Fig [2] shows the NV, when the simulation 
was carried out blocking the destructive effect of noise us- 
ing conditions given in Section HTTl [Eqn. [T] . It shows that 
initially, NV decreases rapidly with noise lavel and after 
reaching the minima it remains constant. As low value of 
NV indicates coherent oscillations, the system remains 
in its coherent state even at higher noise level. This phe- 
nomenon may be termed as constant coherence resonance 
(CCR) as coherence remains unaltered at higher noise 
level. Though CCR has been rarely observed in experi- 
ments, it is one of the important noise induced phenom- 
ena as both CR and SR can be obtained from the basic 
mechanism of CCR. 

When the output of CCR was directly interfered with 
noise, the conventional CR has been obtained in the sense 
that it has unique minimum for optimum noise level. NV 
decreases for small level of noise and then increases with 
increase in noise level as shown by solid line in Fig[21 Min- 
imum of the curve (solid line) represents occurrence of 
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maximum coherence for optimum noise level. At higher 
level of noise NV increases due to the destructive effect 
of noise. For even larger noise the output is totally gov- 
erned by noise and the jumps across the threshold are 
so frequent that the NV decreases again which has also 
been found in noise induced resonances in delayed feed- 
back system [23 |. 



be explained [l^iUHl- In our opinion, adding noise to any 
of the two equations basically destroys the actual dynam- 
ics and hence one does not get CCR. Using our approach 
we have got rid of this problem. 



B. Stochastic resonance 
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FIG. 3: Frequency of the limit cycles with noise: obtained 
from refractory time (dashed line); from interpeak distances 
(dot line); and from FFT (solid line). 

In order to study the frequency change of the au- 
tonomous dynamics in CCR with noise level, we have 
estimated frequency using different techniques. The fre- 
quency obtained from the peak to peak distances, fast 
Fourier transform (FFT) and refractory time with noise 
has been shown in FigOby using dotted, solid and dashed 
lines respectively. Though at small noise level, frequency 
obtained using above three techniques shows some dif- 
ference, at higher noise level they are almost identical. 
From Figs [T] and [31 it is clear that the trend in increase 
in frequency due to noise and due to increase in the CP 
are almost identical. This may be due to increase in ef- 
fective value of CP by noise lavel. 

The experimental results consistent with the above nu- 
merical results have been observed in the unijunction 
transistor relaxation oscillator (UJT-RO) and plasma [1, 
i, il. In UJT-RO, CCR appeared due to the blocking 
of noise from its destructive effect by UJT itself. De- 
tail of the experiment will be found in Ref [3| . Whereas, 
in plasma, sheath was found to be a natural blocker of 
noise [1, In these experiments, the frequency of the 
oscillations was also observed to increase with the noise. 
As the frequency of the autonomous dynamics in these 
cases increased with increase in CP, one can assume that 
the effective increase in CP was due to noise and this is 
consistent with numerical simulation. 

Earlier, in order to show CR in FitzHugh-Nagumo 
model, noise was added to the second Equation [of the 
Eqns [5] which contains the control parameter and role of 
noise was interpreted as an irregular modulation of the 
bifurcation parameter a that switches the limit cycle on 
and off. But adding noise to the first equation also pro- 
duces CR [Eqns U , where switching mechanism cannot 



In order to get stochastic resonance (SR) in the Fitz 
Hugh-Nagumo model, the CP has been perturbed by 
both a subthreshold periodic signal and noise. The main 
mechanism of the SR in an excitable system is that when 
a subthreshold signal is applied along with noise, system 
crosses the threshold at a place of occurrence of the peak 
of the subthreshold signal more frequently and hence sys- 
tem traverses one limit cycles at this place, i.e., the sys- 
tem output mimics the subthreshold signal [1, H^. Here, 
the important point is that the structure of the limit cy- 
cles is determined by the system dynamics and the time 
intervals between them are determined by interpeak dis- 
tances of the applied subthreshold signal. To study SR, 
'a' has been set to oq = 1.25 so that it shows fixed point 
behavior. A subthreshold pulse of the frequency, dura- 
tion and amplitude (A) 0.1, 0.05 and -0.21 respectively, 
is applied along with noise ^ of standard deviation D. 
Whenever, a = A + Di^ + ao < 1 is achieved, the system 
will traverse one limit cycle before coming to its fixed 
point state. The simulation is identical as for the CCR as 
described above except that the CP is perturbed with a 
subthreshold signal and noise. For low level of noise, cor- 
respondence between output and the subthreshold signal 
is very little, and for optimum noise level the correspon- 
dence is excellent and at higher level of noise the system 
response is dominated by noise. 
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FIG. 4: (a) Cross correlation (C) vs D and (b) AMD vs D 
estimated from the time series for stochastic resonance. 

The SR has been quantified by two stochastic param- 
eters absolute mean difference (AMD) 3] and correla- 
tion coefficient(C) !2Qj. AMD is defined as AMD = 
abslmeanitp/S — 1)], where, S is the mean interpeak dis- 
tance of the applied subthreshold signal. C is defined as 
C = mean[{x — mean{x)){y — mean(y))], where, x and y 
are the subthreshold and output signal respectively. For 
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SR, AMD and C should have minimum and maximum 
respectively. Fig 21^ a) and (b) shows C and AMD for an 
applied subthreshold signal of frequency 0.1. Minimum 
in AMD with noise and maximum in C is a typical signa- 
ture of SR. The minimum corresponds to optimum noise 
level at which output mimics the subthreshold signal in 
the maximal manner. The experimental results consis- 
tent with the above numerical analysis have been found 
in the plasma and electro chemical system and details 
will be found in Ref. 

V. CONCLUSION 

Here, we have presented an alternative modeling tech- 
niques to explain a new noise induced phenomenon CCR 
in an excitable system. This also explains other two noise 
induced resonance phenomena CR and SR. The excita- 
tion mechanism which was not clear in the earlier ap- 



proach is now very much clear. Though SR and CR have 
been observed in many excitable system, there are a very 
few experiments on CCR. Possible reason for not get- 
ting CCR in the earlier experiments was that the system 
configurations were unable to block the destructive effect 
of noise. We hope that the experimental observation of 
CCR is possible in the earlier experiments, if the noise is 
allowed such that it only perturbs control parameter by 
choosing suitable parameters and system configurations. 
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